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ABSTRACT 

We investigate gas inflow and outflow histories in Milky Way-like disk galaxies, to get new insights 
into the baryonic processes in galaxy formation and evolution. For this purpose, we solve the equations 
for the evolutions of the surface mass densities of gas and metals at each radius in a galactic disk, based 
on the observed structural properties of distant star-forming galaxies, including the redsliift evolution 
of their stellar mass distribution, their scaling relation between the mass of baryonic components, 
star formation rate (SFR) and chemical abundance, as well as the supposed evolution of their radial 
metallicity gradients (RMGs). We find that the efficiency of gas inflow for a given SFR decreases with 
time and that the inflow rate is always nearly proportional to the SFR. For gas outflow, although its 
efficiency for a given SFR is a decreasing function of time, similarly to gas inflow, the outflow rate is 
not necessarily proportional to the SFR and the relation between the outflow rate and SFR strongly 
depends on the evolution of the adopted RMG. We also find that the results on the outflow rate can 
be reproduced in the framework of momentum-driven (energy-driven) wind mechanism if the RMG is 
steepening (flattening) with time. Therefore if the well measured RMGs and their evolution for Milky 
Way-like galaxies are obtained from future observations, then our results will be useful to constrain 
the main driving mechanism for their galactic outflows. 

Subject headings: galaxies: abundances - galaxies: evolution - galaxies: formation 


1. INTRODUCTION 

Galactic inflow and outflow processes provide a signifi¬ 
cant influence on the budget of baryon and star formation 
activities in galaxies and thus play an important role in 
their formation and evolution. However, how baryonic 
matter is actually funneled into or ejected from star¬ 
forming regions in disk galaxies, like the Milky Way, is 
not clearly understood. 

In a grand picture of galaxy formation and evolution 
provided by ACDM cosmology, baryonic inflow into a 
galactic disk occurs when gas accretes from the inter- 
galactic medium into the dark matter halo and subse¬ 
quently collapses toward its central region. However, 
past semi-analytic models, in which most of the gas in 
the halo are assumed to accrete into the galactic disk, 
predict that the slope at the faint end of the galaxy lu¬ 
minosity function is much steeper than that actually ob¬ 
served (e.g., White & Frenk 1991). This implies that 
the gas cooling and accretion into the disk ought to be 
suppressed in real galaxies by feedback processes asso¬ 
ciated with, e.g., UV radiation, supernovae explosions 
and active galactic nuclei. These feedback processes are 
also necessary to reproduce the structural properties of 
disk galaxies. Indeed, earlier simulations with only mod¬ 
est feedback effects suffer from the so-called angular mo¬ 
mentum problem that the simulated disk galaxies have 
too low angular momentum, thereby being more compact 
and denser than the observed ones (e.g., Navarro & Stein- 
metz 2000). Recent high-resolution hydrodynamical sim¬ 
ulations taking into account stronger feedback effects 
based on more realistic physical background have par¬ 
tially resolved such problems in galaxy formation (e.g., 
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Okamoto et al. 2014; Vogelsberger et al. 2014; Sawala 
et al. 2014), but have not been completely successful yet 
(e.g., Sparre et al. 2014; Oman et al. 2015). 

Feedback processes in galaxies are actually important 
as a driver of galactic outflow, which has been observed 
ubiquitously in starburst galaxies at both high and low 
redsliifts (e.g., Shapley et al. 2003; Rupke et al. 2005; 
Tremonti et al. 2007; Weiner et al. 2009). It is especially 
remarkable from the observations of such galaxies that a 
large amount of metals are expelled from their galactic 
disks (e.g., Bouche et al. 2005, 2006, 2007; Zahid et 
al. 2012; Peeples et al. 2014), thereby implying the 
importance of galactic outflow in the chemical evolution 
of both galaxies and inter-galactic medium. 

To get more detailed insights into galactic inflow and 
outflow processes, it is worth investigating the global 
scaling relations of these star-forming galaxies, such as 
their stellar mass - metallicity (MZ) and stellar mass - 
star formation relations. For example, the MZ relation 
shows that at each redshift there is an empirical upper 
limit in metallicity for massive galaxies. This suggests 
that the chemical evolution in these galaxies accompanies 
gas inflow and outflow processes, unless otherwise the 
galactic metallicities, as derived in the closed-box mod¬ 
els, are predicted to increase monotonously with increas¬ 
ing the stellar to gas mass ratio without any upper limits. 
Recently various studies based on numerical simulations 
and/or semi-analytic models have considered these global 
scaling relations to investigate the basic properties of gas 
inflow and outflow in normal star-forming galaxies (e.g., 
Dave et al. 2011; Lilly et al. 2013; Kudritzki et al. 2015; 
Belfiore et al. 2015). 

We note that most of these previous studies set spe¬ 
cific assumptions on the forms of gas inflow and outflow 
rates without explicit reasoning. For example, one of 
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the popular assumptions is that inflow and outflow rates 
are always proportional to star formation rate (hereafter 
SFR) (e.g., Dayal et al. 2013; Yabe et al. 2015), which 
has been adopted in many studies without concrete as¬ 
sessment from the observational data. Therefore, in this 
paper, we investigate the general properties of gas in¬ 
flow and outflow rates, in particular on their relation to 
SFR without prior assumptions. For this purpose, we fo¬ 
cus on Milky Way-like galaxies (hereafter MWLGs), for 
which star formation histories are well investigated by 
the observational studies of distant galaxies. Provided 
that MWLGs follow the observed global scaling rela¬ 
tions in their metallicities at all their evolutional stages, 
we obtain the inflow and outflow histories of MWLGs 
constrained by the observations of distant star-forming 
galaxies and derive their relation to SFR. 

In addition to these global relations provided by spa¬ 
tially unresolved observations, we also utilize the spa¬ 
tially resolved distributions of metals inside distant disk 
galaxies, which have been recently revealed by the ob¬ 
servations with integrated field units (e.g., Cresci et al. 
2010; Troncoso et al. 2014). Thus, combined with these 
data for internal metallicity distributions and their red- 
shift dependence in distant disk galaxies, we derive the 
time evolution of gas inflow and outflow rates and their 
relation to SFR at each radius in MWLGs. 

This paper is organized as follows. In Section 2, we 
present our method to derive the inflow and outflow rate 
densities from the observational results of distant star¬ 
forming galaxies. In Section 3, we show the inflow and 
outflow histories in MWLGs revealed from our analy¬ 
sis. In Section 4, we discuss the physical interpretations 
of the properties of the derived inflow and outflow rate 
densities. Finally, our conclusions are drawn in Section 
5. 

Throughout this paper, we use the following cosmolog¬ 
ical parameters: Ho = 70 km s _1 Mpc -1 , Ha = 0.7, and 
H m = 0.3, for all the relevant parts of the analysis. 

2. METHOD 
2.1. Basic Equations 

To investigate the surface densities of gas inflow and 
outflow rates in MWLGs at any time, t , and radius, 
7?, along the galactic disk, denoted as Si n (f, R ) and 
Scut {t, R), respectively, we solve the following equations 
(e.g., Pagel 1989; Chiappini et al. 1997), 

BTi 

-£- = “(I - ^)£sFR + - £ out , (1) 

Ot 

d{Z 2 eas) = YY SFn -Z(l-lZ)Y SFR +Z in Y in -Z out Z out . 
at 

( 2 ) 

Equations © and (ED represent the evolutions of a sur¬ 
face gas mass density, E gas , and of a mass fraction of 
heavy elements, Z , respectively. The first term on the 
right hand side of equation © describes the net gas con¬ 
sumption by star formation, where E§fr and 1Z are the 
surface density of SFR and the mass fraction returned 
back to interstellar medium (ISM) via stellar mass loss, 
respectively, and the description of this term is based 
on an instantaneous recycling approximation. In this 
paper we set 7 Z = 0.45 corresponding to the Chabrier 
initial mass function (Leitner & Kravtsov 2011). This 


return fraction includes the contribution of stellar mass 
loss from low mass stars whose lifetime is as long as the 
age of the Universe. This may be somewhat problematic 
when studying an evolution of galaxy. However, accord¬ 
ing to Leitner & Kravtsov (2011), the return fraction 
derived from stars whose lifetime is shorter than 1 Gyr 
is ~ 0.38, so that the mass returned back to ISM from 
low mass stars, for which an instantaneous recycling ap¬ 
proximation may be invalid, is not significant. Therefore 
we expect that even if we adopt such a delayed stellar 
mass loss in our analysis, our conclusion in this paper re¬ 
mains basically unchanged. The second and third terms 
on the right hand side of equation © denote the con¬ 
tributions of inflow and outflow, respectively. The first 
term on the right hand side of equation ED describes the 
supply of heavy elements newly synthesized in massive 
stars, where Y is the nucleosynthetic yield, meaning the 
mass of heavy elements added into ISM per unit SFR. 
We assume that Y is constant. In this paper we focus 
on the chemical abundance for oxygen, which has been 
estimated in many distant galaxies, and adopt the nu¬ 
cleosynthetic yield of Y = 0.015 (Peeples et al. 2014, 
and references therein). The second term on the right 
hand side of equation ED is the mass of heavy elements 
finally locked up in stars. The third and fourth terms on 
the right hand side of equation @ represent the mass 
injection and ejection of heavy elements associated with 
inflow and outflow, respectively, where Z- m and Z out are 
mass fractions of heavy elements in inflowing and out¬ 
flowing gas, respectively. In these equations, we assume 
that inflowing gas contains no heavy elements and the 
metallicity of outflowing gas corresponds to that of the 
ISM, implying Z ln = 0 and Z out = Z , similarly to most 
of the previous chemical evolution models (e.g., Dayal et 
al. 2013; Belfiore et al. 2015). 

Using equations © and (Oil . we can calculate Yj n (t, 77) 
and E out (f, R) if the three quantities, E gas (f, 7?), 
Esfr (t, R) and Z{t, 7?), are given. To derive these three 
quantities in MWLGs, we adopt the recent observational 
results for distant star-forming galaxies as detailed in the 
next subsections. 

These calculations are carried out in a radial range 
from R = 0 to 10 kpc with a radial grid of AT? = 0.1 
kpc. This outer limiting radius roughly corresponds to 
the apertures of the high-redshift galaxies in their ob¬ 
servations of SFR and chemical abundance (Mannucci et 
al. 2010, and references therein). Thus these observed 
quantities derived within the apertures of galaxies will 
be compared with those we derive within 7? = 10 kpc 
in the model galaxy. For calculating the time evolution 
of the system, we adopt the time grid of At = 100 Myr 
and confine ourselves to the redshift range of 0 = 0 to 2, 
within which all the required observational information 
on the structural and chemical evolutions of MWLGs are 
available. 

We note that equations © and @ do not contain the 
effect of radial redistributions of heavy elements inside 
the galactic disk for the sake of simplicity in our follow¬ 
ing analysis. Such radial redistributions of metals can 
actually occur via radial gas flow associated with gravita¬ 
tional disk instabilities (e.g., Athanassoula 1992; Noguchi 
1998) or re-accretion of metal-enriched gas ejected from 
the galactic disk as outflow (e.g., Oppenheimer & Dave 
2008). We will discuss how these effects would modify 
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our results in Section 4.4. 

2 .2. The Structural Evolution of MWLGs 

To estimate E gas (f, R) and EsFR(t, R) for MWLGs, 
we follow the results of van Dokkum et al. (2013) (here¬ 
after vD13), who investigated the stellar mass distribu¬ 
tion of MWLG progenitors selected based on the abun¬ 
dance matching method from 2 ~ 0 to 2.5. We adopt 
here the average evolution of stellar mass, M*, Sersic in¬ 
dex, n, and effective radius, R e g, of MWLGs obtained 
from the fittings to the results of vD13, 

log(M*/M 0 ) = 10.7 - 0.045^ - 0.13z 2 , (3) 

n = 2.74 -0.1472 - 0.1372 2 , (4) 

R eS = 3.42 - 0.1822 -O.I 82 2 kpc. (5) 

Equation ([3]) implies that MWLGs had already obtained 
80 % of their present stellar mass by 2 ~ 1. A similar 
mass evolution is expected from the analyses of disk stars 
in our Galaxy (Snaith et al. 2014). Equations 0 and 
(0 include both the bulge and disk components without 
decoupling each, so the Sersic index at 2 = 0 is 2.74, 
which is greater than the typical value of a galactic disk 
with n ~ 1. The time evolution of a stellar mass pro¬ 
file led from equations Q and 0 is such that prior to 
2 ~ 1 , the growth of stellar mass proceeds almost simul¬ 
taneously at each radius within the disk, whereas, after 
2 ~ 1 , the growth of stellar mass in the central region 
of the disk gradually decreases with time to make the 
stellar mass profile evolve in an inside-out manner. 

We here assume that all the stellar mass at each radius 
originates only from star formation, whereby there is no 
radial migration of stars in the galactic disk. This as¬ 
sumption allows us to derive Esfr^j -R) simply from the 
stellar mass growth. We then translate Esfr(£, R) into 
E gas (f, R) based on the Kennicutt-Schmidt (KS) relation 
between SFR and surface gas density observed in local 
star-forming galaxies, i.e., Esfr, oc E“ as with a ~ 1.4 
(Schmidt 1959; Kennicutt 1998). This assumption for 
star formation law may be reasonable because various 
studies have confirmed that the KS relation is basically 
satisfied up to 2 ~ 2 (e.g., Genzel et al. 2010). We note 
here the recently updated knowledge that Esfr cannot 
be simply represented by a single power law of E gas but 
correlates rather with molecular gas density more tightly 
than atomic one (Bigiel et al. 2008). Gas in a galactic 
disk actually consists of various phases, such as neutral, 
ionized, molecular and dust components, and each state 
can be associated with star formation in a different way. 
However in this work, we simply define the gas density as 
basically the sum of atomic and molecular gas densities 
to avoid the complexity of the transformation between 
these different states. 

With the above procedure, we are thus able to obtain 
both E gas (f, R) and EsFR(t, R) of MWLGs from the re¬ 
sults of vD13. Recently, based on the similar method to 
vD13, Morishita et al. (2015) also investigated the stel¬ 
lar mass distribution of MWLG progenitors from 2 ~ 0.5 
to 3, but adopting a different stacking analysis, which 
takes into account an position angle and axis ratio of 
each MWLG progenitor. Their results partly disagree 
with vD13 in the point that the mass in the central re¬ 
gion of the MWLG progenitor increases with time even 


after 2 ~ 1. We use their results to calculate E gas (f, R) 
and Esfr(L R) and find that the derived histories of gas 
inflow and outflow of MWLGs do not significantly de¬ 
pend on the difference in the adopted evolution of the 
stellar mass profile. In what follows, we present only the 
results derived from the results of vD13. 

We note here that the assumption of no radial redis¬ 
tribution of stars in the galactic disk may be too simplis¬ 
tic, because various theoretical and observational studies 
have suggested the importance of radial migration of disk 
stars (e.g., Sellwood & Binney 2002; Kordopatis et al. 
2015; Morishita et al. 2015; Hayden et al. 2015). How¬ 
ever, how actually the evolution of a stellar mass profile 
depends on the radial migration process is not clearly 
understood yet. Some specific discussion on the effect 
of radial migration in our analysis may be possible by 
adopting chemo-dynamical models such as in Schonrich 
& Binney (2009) and Kubryk et al. (2014). However 
quantifying the effect firmly is yet difficult because the 
predicted properties of radial migration are significantly 
different between models. Although the subject is im¬ 
portant and worth investigating, we attempt to keep the 
clarity of the following analysis by neglecting this radial 
migration effect and thus avoiding the associated uncer¬ 
tainties. 

2.3. The Chemical Evolution of MWLGs 

How the chemical evolution in MWLGs proceeds can 
be inferred from the metallicity measurements for distant 
star-forming galaxies. Mannucci et al. (2010) showed 
that at 2 < 2.5 the average chemical abundances of such 
galaxies can be described as the following function of 
their stellar mass and SFR, 

12 + log(0/H) = 8.90 +0.39a;-0.20:r 2 

( 6 ) 

- 0.077a; 3 + 0.064x 4 , 

where x = log(M*/M 0 ) — O.321og(SFR/M 0 yr -1 ) — 10. 
This tight relation, called the fundamental metallicity re¬ 
lation (hereafter FMR), implies that most of star-forming 
galaxies follows the FMR at all epochs since 2 ~ 2.5. In 
this work, we regard the metallicity derived from the 
FMR as the average metallicity of gas within the galaxy. 
However, it is worth noting that this may not necessar¬ 
ily be a proper approximation, because the metallicity 
based on an analysis of line emissions preferentially re¬ 
flects metallicities of gas in star forming regions, where 
most of photons of useful emission lines in the analysis, 
originate. The relation between metallicity derived from 
integrated spectrum within aperture and that properly 
averaged in the galaxy will be an important subject in 
the future studies with both spatially resolved and unre¬ 
solved data of galaxies. 

Combining the FMR with the stellar mass evolu¬ 
tion of vD13 can provide the total chemical abundance 
and its time evolution in MWLGs, Z tot {t). We de¬ 
fine Z tot = /* max ZT, gas RdR/ / 0 flmax E gas J?di? with i? max 
= 10 kpc, which suggests that Z(t, R) can be calcu¬ 
lated from the evolution of the radial metallicity gradi¬ 
ent (hereafter RMG), dZ(t, R)/dR, along the disk. The 
RMGs of distant star-forming galaxies have been derived 
by dedicated observations with integral field units such 
as VLT/SINFONI and VLT/KMOS. Stott et al. (2014) 
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Fig. 1. — The redshift evolutions of the radial metallicity gradi¬ 
ents for the sRMG (red line) and fRMG (blue line) models. 

showed that the RMGs measured for star-forming galax¬ 
ies at z < 1.5 can be described as an increasing function 
with increasing specific star formation rate (sSFR) as 
follows, 

A1 °^ /H) = 0.023 log(sSFR/yr -1 ) (?) 

+ 0.20 dex kpc -1 , 

namely the slope of the radial metallicity distribution be¬ 
comes steeper with decreasing sSFR. This suggests that if 
MWLGs follow this relation, their RMG gradually steep¬ 
ens with time because sSFR is generally higher at higher 
redshifts. Such evolution of the RMG is also shown in 
Gibson et al. (2013), who simulated the chemical evo¬ 
lution of disk galaxies with strong feedback. Also, the 
RMGs of older disk stars in our Galaxy are found to be 
flatter or positive (Casagrande et al. 2011; Toyouchi & 
Chiba 2014). 

The red line in Figure |T| shows the evolution of the 
RMG derived by combining the RMG vs. sSFR rela¬ 
tion in Stott et al. (2014) with the stellar mass evo¬ 
lution in equation ©. The derived RMG at z = 0 
is —0.053 dex kpc -1 , which is slightly flat compared 
to the present RMG in our Galaxy of A[Fe/H]/AR = 
—0.062 dex kpc -1 (Luck & Lambert 2011). We adopt 
this steepening RMG (hereafter sRMG) with time as a 
standard case to derive Z(t, R), with which we investi¬ 
gate gas inflow and outflow histories in MWLGs. 

Although this sRMG model is actually inferred from 
observations, we also consider the different time evolu¬ 
tion of the RMG in a manner that it is more flattened 
with time, which may be the case if disk galaxies form 
in an inside-out process. This model is referred to as 
the flattening RMG (hereafter fRMG) model in compar¬ 
ison with the sRMG model. To construct this model, 
we make use of the observational fact that the RMGs 
observed in local star-forming galaxies are proportional 
to R-g (e.g., Vila-Costas & Edmunds 1992; Sanchez et 
al. 2014), although the observations of high-z galaxies 
show an opposite result (Stott et al. 2014). If we adopt 
this relation between the RMG and effective radius, the 
RMG gradually flattens with time because the effective 
radius of MWLGs is an increasing function of time. 

The blue line in Figure [T] represents the fRMG 
model derived by combining the relation, RMG = 
RMGq (R e ff {z)/R e ff(z = 0)) -1 , and equation ©, where 



Fig. 2.— The inflow and outflow rate densities (top and bottom 
panels, respectively) as a function of the distance from the galac¬ 
tic center at redshifts of z = 0 (red), 0.5 (green), 1 (blue) and 2 
(purple). The results for the sRMG and fRMG models are shown 
with solid and dashed lines, respectively. 

RMGo is the RMG at z = 0, and is set as the same 
as that for the sRMG model. The comparison between 
these two models for the time evolution of the RMGs al¬ 
lows us to investigate how the difference in it affects the 
estimate of inflow and outflow rates in MWLGs. 

3. RESULTS 
3.1. sRMG Model 

We obtain the inflow and outflow rate densities, E; n 
and Eout, by substituting E gas , Esfr and Z derived in 
Section 2.2 and 2.3 into equations © and ©. Here we 
present the results for the sRMG model. Solid lines in 
Figure [2] show Ei n (top panel) and E out (bottom panel) 
as a function of R at redshifts of 0, 0.5, 1.0 and 2.0. It 
is clear that both Ej n and E out are higher in the inner 
disk region at all epochs and that their central values 
monotonically decline with decreasing redshifts. We also 
find that the radial profiles of E; n and E out evolve roughly 
in an inside-out manner rather than a simultaneous one, 
in contrast to the evolution of the stellar mass profile 
introduced in Section 2.2. 

To get further insights into the properties of gas inflow 
and outflow in MWLGs, we here investigate the depen¬ 
dences of E; n (t, R) and E out (i, R) on Esfr(L R)- This 
procedure allows us to inspect whether the assumption 
adopted in many previous studies, namely that inflow 
and outflow rates in galaxies are proportional to their 
SFR, is indeed valid or not. Solid lines in Figure [3] show 
E; n and E out for the sRMG model as a function of Esfr 
at redshifts of 0, 0.5, 1.0 and 2.0. Black dotted lines 
show Ei n = Esfr (top panel) and E out = Esfr (bottom 
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Fig. 3. — Sj n vs. SgpR (top panel) and S ou t vs. SgpR (bottom 
panel) at redshifts of z = 0 (red), 0.5 (green), 1 (blue) and 2 
(purple). The results for the sRMG and fRMG models are shown 
with solid and dashed lines, respectively. 


panel) for comparison. It follows that the inflow rate 
density roughly satisfies a proportionality to the SFR 
density at all epochs since z = 2. For outflow, while E out 
is roughly proportional to Esfr at high z, this propor¬ 
tionality appears broken at low z. 

To describe these redshift evolutions of the relations 
between the inflow/outflow and SFR densities in a more 
explicit form, we fit E; n and E out at each epoch to the 
following power-law functions of Esfr, 


Si n — ?7inE 


kin 

SFR ’ 


^7out E gpj) • 


( 8 ) 


These parameters, ?^ n /out and fc; n / out , correspond to the 
intercepts and slopes of the profiles shown in Figure [31 
respectively. In the top panel of Figure [4[ we show the 
redshift evolutions of rji n and r) out with red and blue solid 
lines, respectively, indicating that rj i n and ?y 0 ut increase 
with increasing redshift. Also, rj- m is always larger than 
?7outj and is larger than unity at z > 1.5, so that the 
inflow rate dominates the SFR at this early epoch. The 
bottom panel in Figure [4] shows the redshift evolution of 
fci n and fcout with red and blue solid lines, respectively. 
We find that k- ln decreases with time, but is nearly unity 
at all epochs, implying the proportional relation between 
the inflow rate and SFR as mentioned above. In contrast 
to the evolution of k ln , fc ou t evolves from ~ 1.0 at z = 2 to 
~ 0.5 at z = 0. This result suggests that the proportion¬ 
ality of the outflow rate to the SFR as adopted in some 
of previous galactic evolution models is not necessarily 
valid. 


3.2. fRMG Model 


Fig. 4. — The redshift evolutions of the parameters, rj in / out (top 
panel) and fc; n / ou t (bottom panel), defined in equation Isll . The 
results for gas inflow and outflow are shown with red and blue 
lines, respectively, where the solid and dashed lines correspond to 
the results for the sRMG and fRMG models, respectively. 


The results for this model are shown with dashed lines 
in Figures EJ [3] and [4] in comparison to the sRMG model 
with solid lines. We find that although the general evo¬ 
lution of E; n and E out for the fRMG model are similar to 
those for the sRMG model, there is a remarkable differ¬ 
ence between the two models that the radial dependence 
of Ej n and E out for the fRMG model shows no signifi¬ 
cant evolution and remains simultaneous at each epoch, 
in contrast to the inside-out evolution of Ej n and E out 
for the sRMG model. 

These common or different properties of Ej n and E out 
between the two models can be understood in terms of 
the redshift evolution of ?7in/out and &m/out- The top 
panel of Figure [4] shows that for the fRMG model both 
r/i n and y out decrease with time and r/i n is always larger 
than r/ out . These properties of r/i n and r/ ou t for this model 
are the same as those for the sRMG model, implying 
that such properties may reflect the evolution of the to¬ 
tal stellar mass or entire chemical abundance in MWLGs 
rather than the evolution of the RMG in the galactic 
disk. In fact, a similar result was reported by Yabe et al. 
(2015), who investigated the general properties of galac¬ 
tic inflow and outflow by analyzing the scaling relations 
of star-forming galaxies between stellar mass, gas mass 
fraction and chemical abundance at each redshift. How¬ 
ever, there is a difference between the results of the two 
models that r ); n and r ] out for the fRMG model are much 
smaller than those for the sRMG model at high z. This 
may be because for the fRMG model the suppression ef¬ 
fect of chemical evolution due to gas inflow and outflow 
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is smaller than that for the sRMG model to make the 
chemical abundance of the central regions of the disk re¬ 
main higher at earlier epochs, so that the RMG is made 
flatter at subsequent epochs. 

From the bottom panel in FigurelH we can confirm the 
proportionality of the inflow rate to SFR in the fRMG 
model, i.e., Ei n oc Eg‘j; R with fc; n ~ 1 , as also seen in the 
sRMG model. The outflow rate is roughly proportional 
to the square root of the SFR at all redshifts, E out oc 
£gp u R with k ou t ~ 1/2, in contrast to the sRMG model. 
This difference in the redshift evolution of fc out between 
two models reflects the difference in the adopted RMG 
evolution, which will be discussed in detail in Section 4.1. 

4. DISCUSSION 

As presented in the previous sections, we parameterize 
the gas inflow and outflow rates in MWLGs in terms 
of rjin/out and ^in/out as described in equation ([5]). The 
properties of these parameters are summarized as follows. 

• Both 77i n and 7f out are a decreasing function of time. 

• The index fc; n is almost unity at all epochs, im¬ 
plying that the inflow rate can be regarded to be 
proportional to the SFR. 

• The index fc out for the fRMG model is about 0.5 
at all epochs, whereas that for the sRMG model 
evolves from ~ 1 at z = 2 to ~ 0.5 at z = 0. There¬ 
fore the evolution of fc ou t is expected to depend on 
the adopted evolution of the RMG. 

These results reflect basic baryonic physics associated 
with gas inflow and outflow processes in the disk. In this 
section, we discuss how the above properties of ? 7 in/out 
and fci n / ou t are related to and deduced from such baryonic 
physics. 

4.1. The Relation between the Evolutions of k ou t and 
the RMG 

The most remarkable difference between the results of 
the sRMG and fRMG models is the redshift evolution of 
kout ■ This difference suggests that the dependence of out¬ 
flow process on star formation activity affects the chem¬ 
ical evolution over the galactic disk significantly. The 
relation between fc out and the RMG can be simply inter¬ 
preted as follows. 

For an ordinary galactic disk that has more active star 
formation in its inner region, fc out > 1 means that the ra¬ 
tio of the outflow rate to SFR becomes higher toward the 
inner region. Therefore, the enrichment process driven 
by active star formation in the inner region of the disk 
is suppressed by the effect of much higher outflow rate 
of gas, so that the steepening of the RMG with time 
is suppressed, although the actual situation is expected 
to depend on the nucleosynthetic yield and initial dis¬ 
tribution of heavy elements in the disk. In contrast, for 
fcout < 1 the chemical evolution always progresses faster 
in the inner disk region, so that the RMG is steepened 
with time. Therefore, for the sRMG model fc ou t > 1 at 
high z leads to the flat RMG and the gradual decrease 
of fc ou t accompanies the steepening of the RMG. On the 
other hand, for the fRMG model fc ou t -C 1 at all redshifts 
makes the RMG steeper at earlier epochs. 


Thus the chemical evolution over the galactic disk 
strongly depends on the properties of the outflow pro¬ 
cess, whereby the physical origin of the evolution of fc ou t 
is a key to understanding the chemical evolution of the 
galactic disk. In the next subsection, we present possi¬ 
ble models for explaining these properties of the outflow 
process. 

4.2. The Relation between E ou t and Esfr 

Previous studies with both analytical models and nu¬ 
merical simulations have suggested that gas outflows 
from a galactic disk are driven by the injection of ki¬ 
netic energy (e.g., Springel & Hernquist 2003; Okamoto 
et al. 2010) or momentum (e.g., Murray et al. 2005; 
Oppenheimer & Dave 2006; Hopkins et al. 2012) from 
stellar feedback into ISM. The former and latter outflow 
mechanisms are called the energy-driven wind (EDW) 
and momentum-driven wind (MDW), respectively. Re¬ 
cently Okamoto et al. (2014) showed that a baryonic 
feedback model including both EDW and MDW success¬ 
fully reproduces the various properties of galaxies, for 
example stellar mass function, stellar mass - halo mass 
relation and their redshift evolutions from z = 4 to 0. 
However, which wind mechanism is actually the main 
source of galactic outflow is not well understood. In this 
section we consider these EDW and MDW models as the 
origin of the relation between E out and Esfr obtained 
in our analysis, which will provide valuable insights into 
the main mechanism of galactic outflow. 

4.2.1. Energy-driven wind 

First, we attempt to describe E out as a function of 
Esfr in the framework of the EDW model. Here we 
consider a star cluster with stellar mass of M c \. Because 
the energy injected from the star cluster into ISM is pro¬ 
portional to M c i in the EDW model, the mass of outflow 
ejected from the star cluster, m w , is described by using a 
typical velocity of outflowing gas when leaving the disk, 
u w , as follows, 

M cl , 

m w oc . (9) 

V 

We assume that v w corresponds to the escape velocity of 
the disk defined as, 

u w ~ \J TrCTJ(E star + E gas ) oc ^/(l + /r-^Egas , (10) 

where /i s E gas /E sta r, and H is the scale height of the 
disk and is assumed to be constant along the radius. 
As the number of star clusters with M c \ in the region 
where stars are formed in the rate of Esfr is estimated 
as Esfr/A/cU we can describe E out in terms of Esfr as 
follows, 

Eout ~ W w S . S > FR ■ (11) 

M cl 

By substituting equations © and m into m and 
adopting the KS relation, Esfr oc E“ as , we express E out 
as functions of fj, and Esfr as follows, 

V nr ^ y( Q—1 )/ a /I O'! 

^out OC -^SFR • \ LZ ) 

_L r fJj 

We regard the factor /x/(l + /r) in this equation depends 
only on time, because the variation of fi along the radius 
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Fig. 5.— The redshift evolution of r/ ou t divided by ///(1 4- //) 
and /j,/( 1 + //) for the sRMG model (red solid line) and fRMG 
model (blue solid line), respectively. The parameter // used here is 
calculated from the total gas and stellar mass at each redshift. For 
comparison, we also show the undivided jj ou t for the sRMG model 
(red dashed line) and fRMG model (blue dashed line). Each line 
is normalized by the value at z = 0. 

is much smaller than that of Esfr- Therefore /x/( 1 + \x) 
and (ct — 1 )/a in this equation correspond to rj on t and 
fc out in equation ID, respectively. 

This formulation reveals that k ou t is independent of 
time and that by adopting a = 1.4 as observed in local 
star-forming galaxies we find k out = (a — l)/a ~ 0.29. 
These properties of fc out are basically in agreement with 
our results obtained for the fRMG model. Moreover, 
since p is generally much higher at higher redshifts, we 
reproduce ?; 0 ut as a decreasing function of time. To con¬ 
firm the validity of equation m, in Figure [5] we show 
?7out divided by the factor n/(l + p) with blue solid line 
for the fRMG model. In this figure we also show the un¬ 
divided ^out with blue dashed line for comparison. We 
note that each rj out is normalized by the value at z = 0, 
and /x used in this procedure is derived by using the total 
gas and stellar mass at each redshift. We find that the 
divided r/ out for the fRMG model is only slightly smaller 
than unity at all redshifts and is approximately constant 
at z > 0.25. Thus the EDW model can well represent the 
properties of gas outflow obtained for the fRMG model. 
In contrast, the properties of fc out in the sRMG model 
are not understood in the framework of the EDW model 
but the MDW model is at work as explained below. 

4.2.2. Momentum-driven wind 

We now derive the relation between E out and Esfr 
based on the MDW model. We consider that the radia¬ 
tion emitted from a star cluster drives the ambient gas as 
outflow. The momentum injected into the ambient gas 
from a star cluster via the radiation during the time scale 
of 7i n j, which should roughly correspond to the lifetime 
of massive stars, i.e., ~ 10 Myr, is described as follows, 

P ~ (1 + Tir) — T in j , (13) 

where L and tir are the luminosity of the star cluster 
and the optical depth of the ambient gas to the infrared 
emission by dust grains, respectively. 

We note that tir is proportional to the surface gas den¬ 
sity of giant molecular clouds (GMCs) surrounding the 
star cluster, Egmc, rather than E gas as in equations ID 


and ©■ Therefore, in order to discuss the result in Fig¬ 
ure [4j we need to explicitly describe the relation between 
Egmc and E gas . We consider here that each of a GMC 
is originally a fragment of the galactic gas disk with E gas 
and the scale of fragmentation is roughly consistent with 
the scale height of the gas disk. Then the mass of a GMC 
can be described as Mqmc = 7 rE gas lL 2 = 7tEgmcFgmC’ 
where tgmc is the radius of a GMC. Assuming that the 
ratio of H to tgmc is nearly constant, as in Murray et 
al. (2011), we obtain tir oc E gas . In this manner, we 
connect the physics in a galactic disk scale with that in 
a GMC scale in this discussion. 

Here we consider the optically thin (tir <C 1) and thick 
(fir 1 ) limits, and the mass of outflow ejected from 
the star cluster for each limit is described as follows, 

f Finj (T/c)/u w OC M c \/v w 

(optically thin limit) , 
nni(riRL/c)/v w oc E gas M c i/u w 

(optically thick limit) , 

where we assume that L is proportional to M c \. From 
equations (flOl) . (fill) and (fill) . E out for the MDW model 
can be presented as follows; 

(^sfr 1 ^ 2a (optically thin limit) 

S ° Ut * \ \[rh s sf“r 1)/2 “ (optically thick limit/^ 

Thus, we find that E out for the MDW model can also be 
expressed in the same functional form as equation m, 
similarly to the EDW model. 

The remarkable difference from the case of the EDW 
model is that k ou t depends on the optical depth of the 
ambient gas and thus the surface gas density. Adopting 
a = 1.4 leads to fc ou t ~ 0.64 and 1.36 for the optically 
thin and thick limit, respectively, which roughly repro¬ 
duce the values of k out for the sRMG model at z ~ 0 and 
2, respectively. Therefore, we suggest that the evolution 
of fcout obtained for the sRMG model is understood in the 
framework of the MDW model, as the transition from the 
optically thick to thin limit in the galactic gas disk. This 
transition occurs at the gas density of ~ 100 M 0 pc -2 
(Murray et al. 2011), which would be sufficiently lower 
than the gas densities in high -z MWLGs. 

Equation (fl5l) also predicts the decreasing ri out with 
time. The red solid and dashed lines in Figure [D show 
77out for the sRMG model divided and undivided by the 
factor yjp/{ 1 + p), respectively. From this figure, we 
find that the divided ?y 0 ut remains nearly unity and con¬ 
stant in comparison with the undivided xyout- Thus MDW 
appears to explain well the properties of gas outflow in 
the sRMG model. 

From the above discussion, we conclude that clarifying 
how the time evolution of the RMG in MWLGs pro¬ 
ceeds, i.e., more steepened or flattened with time, dis¬ 
tinguishes which wind-blowing mechanism, injection of 
momentum or kinetic energy from stellar feedback, is 
actually at work. Thus, future observational constraints 
on the RMG in MWLGs will be important to understand 
the main mechanism driving galactic outflows. 

4.3. The Proportionality of Ej n to Esfr 
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In Section 3, we have shown the linear relation between 
E; n and Esfr, which is almost independent of the differ¬ 
ence in the RMG evolutions. A similar proportionality 
was also confirmed by Recchi et al. (2008), who com¬ 
pared a chemical evolution model including gas inflow 
being proportional to SFR with a model assuming expo¬ 
nentially decreasing inflow rate with time. We show here 
that this proportionality is indeed universal for all the 
galactic evolution models with various gas inflow histo¬ 
ries. 

Following the results of this work as well as previous 
works, gas inflow is always systematically larger than 
outflow, thereby we assume here E; n — E out ~ E; n as 
a working hypothesis in the following discussion. Then 
equation m for the conservation of gas mass is rewritten 
as, 

E in = ^jp + a-T^EsFR . (16) 

Since dE gas /df is described with /i as follows, 


dE g as v dln/i 

= £gas—+ Ml - Tw^SFR j 

then, by substituting this equation into equation 
E; n can be expressed as, 


(17) 


csD, 


Ein = |(1 + M(1 ~ K) + d | 1 7 / } S S FR ■ (18) 


The second term in the bracket on the right hand side 
corresponds to the ratio of gas consumption timescale, 
Tgas (= Egas/EsFR), to the timescale on which /i is chang¬ 
ing, t m (= dt/dln/a). We note here that depends on 
both SFR and inflow rate, whereas r gas depends only on 
SFR. This implies that while T gas ~ in the closed- 
box model, the model taking into account inflow yields 
Tgas < t /J; , because the gas supply by inflow cancels 
out the decrease of the gas mass due to star formation. 
Therefore we find from equation (fl8l) that Ei n is basically 
proportional to Esfr for fi > 1, which is in particular the 
case at high redshifts. Then the property of rji n is domi¬ 
nated by the factor (1 + /a) (l — 1Z), which is a decreasing 
function of time, thereby reproducing the time evolution 
of ?7in presented in Section 3. For /i < 1, as expected 
at most of the lower redshifts, the second term in the 
bracket on the right hand side of equation (fTSl) provides 
the deviation from the simple proportionality, which may 
reproduce k- m being slightly smaller than unity as derived 
in our analysis. 

Thus the proportionality of Ei n to Esfr as revealed by 
our analysis is simply understood from the conservation 
of gas mass. This proportionality may imply that gas 
inflow regularly controls the gas budget and SFR in the 
galactic disk. 


4.4. The Influence of Radial Redistributions of Metals 
on Our Results 

As noted in Section 2.1, our analysis does not include 
the effect of radial redistributions of metals inside the 
galactic disk. However, the recent observations of nearby 
galaxies having extended gas disks show that the star for¬ 
mation activities observed at their outskirts are too weak 
to produce their metallicities observed in such regions, 
and therefore suggest the need of significant transporta¬ 
tions of heavy elements from their inner to outer disk 


regions (Werk et al. 2010, 2011; Bresolin et al. 2012). 
Such transportations of heavy elements may be driven 
by radial flows of gas inside galactic disks due to grav¬ 
itational instabilities or re-accretion of metal-enriched 
galactic wind into the galactic disk. In this subsection we 
discuss the possible impacts of such radial redistributions 
of heavy elements on our results. 

4.4.1. Radial flows of gas 

Radial flow of gas may be frequently driven by non- 
axisymmetric structures such as bar/spiral and signifi¬ 
cantly affect the structural and chemical evolutions of 
disk galaxies (e.g., Sclionrich & Binney 2009; Mincliev 
et al. 2013; Kubryk et al. 2014). We here attempt to 
simply evaluate to what extent such radial flows can ac¬ 
tually change our results. It is worth noting, however, 
that the answer may sensitively depend on the detailed 
properties of given radial flow. 

The effect of radial flow on the evolution of a sur¬ 
face gas mass density at each radius is represented by 
adding an advection term to the right side of equation 

(1) , —d(RY, gas vn)/RdR, where vr is the velocity of ra¬ 
dial flow. Taking into account E; n !§> E out as obtained 
above, we write 

^P~-(l-ft)E S FR + E in —IA (i? E gas u fl ) . (19) 

Several chemical evolution models for our Galaxy includ¬ 
ing radial flow suggested the advection term in equation 
(HD to be of the order of ~ O.OIM© yr 1 kpc 2 (e.g., 
Schonrich & Binney 2009; Kubryk et al. 2014), which 
is comparable to Ej n at z < 0.5 as followed from Figure 

[2] Therefore the influence of radial flow on the estimate 
of Ei n may be significant only at z < 0.5. We also note 
that in our analysis all the terms except for those of 
E; n and advection in equation (fl9l) are derived indepen¬ 
dently from the latter two terms, so that the change of 
E;„ associated with radial flow simply reflects the sign 
and amplitude of the advection term. 

We now consider the effect of bar/spiral structures, 
which generally drive radial flow with negative (positive) 
vr for the gas at R smaller (larger) than the co-rotation 
radius of bar/spiral, R co . Then the last term of the right 
hand side of equation (flTil) is negative at R around R co , 
thereby to compensate this export of gas due to the radial 
flow, Ein becomes larger than that for the case without 
radial flow. In contrast, the advection term is positive in 
both the central and outer disk regions, so Ei n becomes 
lower. Then k- m in the case of non-zero vr at R < Rco 
(R > R co ) becomes smaller (larger) than that in the case 
of vr = 0. Thus radial flow will provide non-monotonous 
changes in the relation between Ei n and Esfr- However 
when considering the effects of radial flow over the entire 
disk, the difference of k m between the case of zero and 
non-zero vr may be small, implying that Ei„ is on av¬ 
erage proportional to Esfr even including the effect of 
radial flow. Thus the effect of radial flow on the derived 
properties of gas inflow is expected to be insignificant. 

The effect on the outflow rate, E out , by the presence of 
radial flow can be inferred as follows. For the disk with 
the negative RMG, because metal-poor (rich) gas present 
in the outer (inner) disk regions flow inward (outward), 
the RMG tends to be flattened by such radial flow. Then 
fc ou t may become smaller in response to this flattened 
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RMG since the outflow with smaller k out works so as 
to steepen the RMG more efficiently, as mentioned in 
Section 4.1. 

4.4.2. Re-accretion of metal-enriched galactic wind 

Recent hydro-dynamical simulations for galaxy evolu¬ 
tion suggest that some fraction of gas ejected from the 
galactic disk by feedback processes remain in the host 
halo and accrete back into the disk (Oppenheimer & 
Dave 2008; Dave et al. 2011). Moreover outflowing gas 
may acquire additional angular momentum via feedback 
processes and preferentially re-accrete into outer disk re¬ 
gions than the radii where they originally resided (Bekki 
et al. 2009). Therefore re-accretion of metal-enriched 
galactic wind, if it occurs, may play a role in transport¬ 
ing heavy elements from the inner to outer disk regions, 
thereby strongly affecting the radial distribution of heavy 
elements inside the galactic disk (e.g., Tsujimoto et al. 
2010; Gibson et al. 2013). 

To investigate this influence on our results, we consider 
the case that inflowing gas consists of not only primor¬ 
dial one, but also re-accreting enriched one, for which we 
adopt a non-zero Z- m to solve equation ©. For this ex¬ 
periment, we assume that the metallicity of inflowing gas 
at each radius is proportional to the total metallicity of 
the galaxy, Zi n (t, R ) = a(R)Z to t(t ), where the propor¬ 
tional coefficient, a(R), is assumed to increase linearly 
along the radius so as to represent re-accretion of metals 
preferentially into the outer disk regions. As a fiducial 
example, we adopt a(R) = 0.1 + 0.1(i?./10 kpc) so that it 
changes from 0.1 at R = 0 kpc to 0.2 at R = 10 kpc. We 
believe that this case provides sufficiently large Z- m to 
examine the effect, because the actually observed metal- 
licities in high velocity clouds presently falling into our 
Galactic disk are estimated as only 0.09 times the solar 
metallicity (e.g., Wakker et al. 1999). 

We perform this calculation for both the sRMG and 
fRMG models and find that this form of non-zero 
Zi n (f, R) makes about 20 % of metals ejected from the 
entire disk eventually fall back at later epochs. We also 
find that the results from these calculations are basically 
the same as those shown in Section 3: the change in the 
parameters i] i n /out an d fcm/out by adopting this non-zero 
Zi n are typically only about 10 % and generally con¬ 
fined to less than 30 % at all epochs for both the sRMG 
and fRMG models, so the redshift evolutions of these 
parameters remain the same as those shown in Figure [l] 
Therefore we conclude that the effect of re-accretion of 
metal-enriched galactic wind does not significantly mod¬ 
ify our conclusions in this paper. 

4.5. The Total Inflow and Outflow Rates 

The total inflow and outflow rates at each redshift are 
obtained by integrating the inflow and outflow profiles 
along the radius shown in Figure [2] The information 
on these total gas flow rates, which are to be compared 
with the observations of distant MWLGs, can place use¬ 
ful constraints on galactic chemical evolution models. We 
here present the integrated inflow and outflow masses ob¬ 
tained from our analysis, and discuss their roles in galaxy 
evolution. 

4.5.1. Inflow 



Fig. 6. — The total inflow rate obtained by the integration of the 
inflow profile along the radius within R = 10 kpc (solid line) and 
the cumulative inflow mass after z = 2 (dashed line) as a function 
of redshift. The red and blue lines show the results for the sRMG 
and fRMG models, respectively. For comparison, the case for the 
ideal baryonic accretion is shown with black lines. 

Figure [ 6 ] shows the total inflow rates (solid lines) and 
cumulative inflow masses after z = 2 (dashed lines) de¬ 
rived by integrating the inflow profile within R = 10 kpc, 
which is the outer boundary of the disk in our analysis as 
mentioned in Section 2, where the results from the sRMG 
and fRMG models are shown with red and blue lines, re¬ 
spectively. We find that the total inflow rates are highest 
at z = 2 and decrease rapidly with time at later epochs. 
It is clear that the total inflow rates for two models are 
remarkably similar, although the internal properties of 
gas inflow between the models are largely different. This 
result suggests that the total mass of inflowing gas into a 
galaxy is determined by the cumulative properties of the 
galaxy, such as its total stellar mass, SFR and mass bud¬ 
get of heavy elements, rather than the spatially resolved 
properties, such as the radial distributions of baryonic 
components and their metallicities. The total masses ac¬ 
creted on the galactic disk after z = 2 are 5.4 x 1 O 1 O M 0 
and 5.5 x 10 10 Mq for the sRMG and fRMG models, re¬ 
spectively. These masses correspond to the present stel¬ 
lar disk mass in MWLGs. 

In Figure [ 6 ] we also show the ideal baryonic accre¬ 
tion into the halo of MWLGs with black lines, which 
is obtained by multiplying the cosmic baryonic fraction, 
/b ~ 0.17, by the average virial mass growth rate for 
MWLGs (Bouche et al. 2010, and references therein). 
We find that this ideal accretion is clearly different from 
the actual inflows presented by our analysis, as the latter 
decline much faster than the former. This difference im¬ 
plies that not all the gas existing in the halo of MWLGs 
has been supplied on the galactic disk, in particular at 
lower redshifts. This suppression of gas infall may be 
caused by feedback from supernovae or active galactic 
nucleus in the galaxy, which makes halo gas heated up 
before accreting on the galactic disk. Recently, Lu et al. 
(2015) investigated the semi-analytic model, which in¬ 
cludes the effects of preheating mechanisms, and found 
that the model can reproduce remarkably well various ob¬ 
servational scaling relations such as the cold gas fraction 
- virial mass relation, disk size - stellar mass relation and 
its evolution. Thus the feedback processes in a galaxy 
may be linking up not only to the production of outflow, 
but also to the several properties of star-forming galaxies 
including baryonic accretion. 
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Fig. 7.— The total outflow rate obtained by the integration of 
the outflow profile along the radius within R = 10 kpc (solid line) 
and the cumulative outflow mass after z = 2 (dashed line) as a 
function of redshift. The red and blue lines show the results for 
the sRMG and fRMG models, respectively. 

4.5.2. Outflow 

Figure [7] is the same as Figure E3 but for the total 
outflow rates. We find that these rates show nearly the 
same evolution for both models, i.e., declining rapidly 
with time, similarly to the total inflow rates. The total 
masses ejected from the galactic disk after z = 2 are 
3.4 x 1O 1O M 0 and 3.5 x 1O 1O M 0 for the sRMG and fRMG 
models, respectively. 

Figure [5] shows the total outflow rate of metals and 
their cumulative mass. This flow is also decreasing with 
time, but this decreasing rate is somewhat slower than 
that for gas outflow because the chemical abundance of 
the galactic disk is a increasing function of time. The to¬ 
tal masses of metals ejected from the galactic disk after 
z = 2 are 4.7 x 1O 8 M 0 and 4.3 x 1O 8 M 0 for the sRMG and 
fRMG models, respectively. These masses correspond to 
about 40 % of the mass of all metals synthesized in the 
star formation history of MWLGs, Y(1 — 1Z)M *. Similar 
results were reported by the previous works, which com¬ 
pared the mass of metals still contained in the galactic 
disk with that produced by all the star formation pro¬ 
cesses (e.g., Zahid et al. 2012). According to Peeples et 
al. (2014), about 80 % of metals are already lost from 
the galactic disk in MWLGs probably due to outflow. 
Thus these works including our current analysis imply 
that galactic outflows have provided significant effects 
on the chemical evolution of galaxies. 

5. SUMMARY AND CONCLUSIONS 

We have investigated the gas inflow and outflow rate 
densities in Milky Way-like galaxies as a function of the 
SFR density, in order to get new insights into the bary- 
onic processes tightly relating to galaxy formation and 
evolution. For this purpose, we solve equations m and 
m for the evolutions of the surface mass densities of gas 
and metals at each radius in a galactic disk, based on 
the observational studies of distant star-forming galax¬ 
ies, such as the redshift evolution of their stellar mass 
distribution, their fundamental metallicity relation and 
the evolution of their radial metallicity gradients. The 
solutions of these equations allow us to evaluate the effi¬ 
ciency of inflow and outflow rate densities, Ei n and S out , 
and their dependence on the SFR density parameterized 
with r] in/out and ^’m/out as presented in equation (0). 



Fig. 8. — The total metal outflow rate obtained by the integration 
of the metal outflow profile along the radius within R = 10 kpc 
(solid line) and the cumulative metal outflow mass after 2 = 2 
(dashed line) as a function of redshift. The red and blue lines show 
the results for the sRMG and fRMG models, respectively. 


We have found that k\ n is approximately unity at all 
redshifts, representing the proportionality of the inflow 
rate to the SFR, and rj- ln is a decreasing function of time. 
Such properties of gas inflow can be understood based on 
the equation of gas mass evolution. Therefore gas inflow 
dominates the budget of gas and star formation activity 
in a galactic disk. 

For gas outflow, while ?y 0 ut decreases with decreasing 
redshift, similarly to ry; n , fc ou t deviates largely from unity 
at most redshifts, indicating that gas outflow rate is 
not necessarily proportional to SFR. Moreover, we have 
found that the relation between the outflow rate and SFR 
strongly depends on the evolution of the adopted radial 
metallicity gradient. Thus galactic outflow provides the 
significant influence on the distribution of heavy elements 
in disk galaxies. 

The relation between the outflow rate and SFR is 
understood based on the momentum and energy-driven 
wind models. It is found that momentum-driven (energy- 
driven) wind mechanisms can well reproduce the red¬ 
shift evolution of k OVLt based on the steepening (flatten¬ 
ing) radial metallicity gradients with time. However, 
the detailed spatial distribution of metallicity in Milky 
Way-like galaxies has not been constrained well from 
the observations of distant galaxies. Future observations 
armed with the state-of-the-art surface-spectroscopic in¬ 
struments such as VLT/KMOS will be important to un¬ 
veil the detailed spatial distribution of metallicity in dis¬ 
tant galaxies at several redshifts. If the well measured 
radial metallicity gradients and their evolution are avail¬ 
able from such observations, then the method that we 
have presented in this work will be useful to constrain 
the main driving mechanism for their galactic outflows. 
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